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PREFACE 

This  work  was  conducted  in  support  of  NUSC  Project  No.  £63400, 
"Time/Angle  Spreading  of  Bottom  Energy",  Principal  Investigator, 

E.  Jensen  (Code  3332).  The  sponsoring  activity  is  NOPDA  Code  323 
(D.  King)  under  program  element  63704N. 
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A  Relationship  Between  off-Specular  Scattering  and 
Time  Spread  Phenomena 


1.  Introduction 


Acoustic  signals  received  from  near  sea  surface  explosive 
detonations  were  observed  to  spread  in  time  of  arrival  due  to 
bottom  interaction.  The  time  spread  extends  considerably  beyond 
the  pulse  width.  In  this  study,  off  specular  scattering  due  to 
bottom  roughness  which  is  thought  to  be  a  possible  mechanism 
causing  the  spreading,  is  investigated.  The  scattering  theory  of 
Marsh-Kuo-1  is  utilized  in  an  attempt  to  predict  the  observed  time 
spreading . 

2 .  Selected  Geometry 


For  many  experiments  the  source  and  receiver  are 
approximately  at  the  same  distance  K  from  the  bottom.  For  this 
geometry,  off-specular  arrivals  at  the  same  time,  t  ,  are 
reflected  from  an  elliptical  ring: 

<& /!>$■"  1  ,  ( 


where 


a1  =  c+t+I^ 


(2) 


b’  =  ect hftt'-txv  e  , 

"tj  ■=.  specular  arrival  time, 
ia  -  direct  arrival  time, 

C  =  acoustic  speed. 


3 .  Model  Description 

Utilizing  a  method  of  perturbation,  Kuo*  generalized  the 
scattering  theory  of  Marsli  and  predicted  reflected  field 
characteristics  from  a  rough  transmitting  bottom.  To  the  second 
order  of  perturbation,  the  reflected  field  consists  of  specular 
and  scattered  components. 

The  specular  component  is  a  concentrated  beam  in  a  forward 
specular  direction.  Its  energy  level  depends  on  Rayleigh' s 
reflectivity  and  the  amount  of  scattered  energy  lost  in  all  other 
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directions.  The  scattered  energy  lost  in  all  directions  depends 
on  the  entire  range  of  the  bottom  roughness  spectrum. 

The  scattered  component  is  diffused  and  its  energy  level  in 
any  direction  is  dependent  on  a  particular  part  of  the  bottom 
roughness  spectrum.  Therefore,  the  scattering  process  in  an 
off-specular  direction  is  roughness  wave  number  selective. 

Given  locations  of  the  point  source  and  receiver  as  indicated 
in  Figure  1,  each  ray  emerging  from  the  source  is  assumed  to  be 
approximately  plane  at  the  scattering  boundary.  The  first 
arrival  after  the  direct  arrival  is  the  specular  component  which 
is  reflected  at  the  boundary  point  with  coordinates  (0,0,0).  The 
succeeding  arrivals,  at  a  particular  time,  are  all  off-specular 
components  arriving  simultaneously  from  a  particular  elliptical 
ring  depicted  in  the  same  figure.  Figure  1  indicates  one  such 
off-specular  component  arrival.  To  find  the  received  level  at  a 
particular  time,  one  must  sum  up  all  contributions  from  such 
arrivals  along  a  particular  elliptical  ring. 

Kuo's^  scattering  strength,  m" ,  predicts  off-specular 
components  to  be  proportional  to  the  product  of  the  scattering 
function  and  the  roughness  spectral  level  at  a  selected  wave 
number.  The  scattering  function  is  independent  of  the  roughness 
but  is  dependent  on  the  incident/reflected  wave  and  media 
properties.  The  selected  wave  number  is  the  difference  wave 
number  between  incident  wave  number  t$,  =.  (zir  (i  ) 

and  the  reflected  wave  number  =  C2iri /c  )  <  \,>*0 
where  ^  =  acoustic  frequency,  C*,(3>  =  incident  wave 

direction  cosines  with  respect  to  (  x,  ^  ),  (  )=  reflected 

wave  direction  cosines  with  respect  to  (  ).  The  magnitude 

of  the  selected  wave  number  can  be  expressed  in  terms  of 
geometrical  quantities  as  given  by  the  following  expression. 


tirM®  viV*.  (rf  “  t1}’- 

(4) 

+  *  {Vfcf  J*- 

where  y-i.  yf 

(5) 

(6) 

-  ”217^  / C 

x,  y:  satisfies  equation  (1). 

The  selected  wave  number  given  by  (4)  will  vary  from  zero  at 
specular  point  (0,0,0)  to  2k^  at  eo  .  The  latter  selected  wave 
number  is  also  the  selected  wave  number  in  the  back-scattering 
direction.  The  extent  of  variation  in  the  selected  wave  number 
along  each  ellipse  was  investigated.  It  was  found  that  the 
magnitudes  of  the  selected  wave  numbers  along  any  ellipse  vary  as 
much  as  a  factor  of  2  within  the  area  between  the  source  and 
receiver.  However,  outside  of  that  area  variations  were  small. 
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Received  level  (RL)  can  be  expressed  in  terms  of  the 
scattering  strength  level,  10  log  m" ,  by  the  following 
expression : 


RL  *  10  log  m"  +  SL  -  TL  +  10  log  A 


(7) 


,  where  SL  =  source  level, 

TL  =  total  transmission  loss  level, 
A  =  insonified  area. 


For  negligible  water  medium  absorption  loss,  the  total 
transmission  loss  becomes  the  double  spherical  spreading  loss 
given  by  the  following  expression. 


TL  =  20  log  r^  +  20  log  rj 


(8) 


The  received  level  at  a  given  time,  ,  is  considered  to  be 
contributed  to  by  a  sonified  elliptical  ring  area  of  one  pulse 
width.  The  corresponding  selected  wave  number  contained  in  the 
scattering  strength  varies  over  this  insonified  ring  area  and 
the  effect  of  each  elementary  area  on  RL  should  be  integrated 
over  the  entire  elliptical  ring.  However,  for  this  study,  a  much 
simpler  approximate  method  described  below  is  utilized. 

4 .  Results  for  a  Rigid  Boundary 

The  mathematical  model  presented  by  Kuo^  in  1964,  was 
intended  for  scattering  from  a  non-shear  supporting  sedimentary 
layer  of  an  infinite  thickness.  As  such,  the  model  can  not 
really  predict  scattering  from  basement  material.  A  model 
development  for  shear  supporting  sediment  and  basement  layers  is 
underway.  One  approximation  that  can  be  made  at  this  time  is  to 
let  the  density  of  the  sediment  become  very  large  in  comparison 
with  that  of  water.  Effectively  this  causes  the  sediment  to 
approach  a  rigid  boundary  and  simulates  solid  basement  material 
to  some  extent. 

For  this  approximation,  the  scattering  strength  (m" )  is 
simplified  to  the  following  expression: 


.«  2^/  3r3ir\  E 

*  T  V  )  U,-*r| 
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where  E  represents  a  boundary  roughness  spectrum.  In  deriving 
the  above  expression,  bottom  roughness  is  assumed  locally 
isotropic  n^ar  the  selected  wave  number  region  -  on  the  order  of 
0.1  (meter  ).  Spof ford/Holmes/SpiessJ  show  measurements 
indicating  approximate  isotropic  roughness  down  to  wave  number  of 
about  10  (meter  A ) .  For  the  following  calculation,  the  power 
law  spectrugi, 

E(k)  =  EQk  ,m  =  0  to  3,  will  be  assumed. 

Because  of  geometrical  symmetry  and  symmetry  of  ,  and 
in  equation  (9),  scattering  strength  varies  only  in  one  quadrant. 
Scattering  strengths  at  three  locations  in  one  quadrant  are 
calculated  and  averaged.  The  averaged  value  is  then  used  in 
equation  (7)  as  applicable  to  the  entire  elliptical  ring  area. 
Source  level  is  assumed  zero  dB. 

Figure  2  depicts  the  time  spread  in  received  levels 
(propagation  loss)  given  by  off-specular  scattering.  Signal 
frequency  and  pulse  length  are  500  B  and  20  ms,  respectively. 

The  specular  arrival  time,  is  calculated  to  be  13.849  seconds. 

The  first  off-specular  arrival  level  is  calculated  at  30  ms  after 
the  specular  arrival.  The  succeeding  off-specular  levels  are 
calculated  at  20  ms  intervals.  Each  curve  represents  the  result 
of  using  a  different  power  law  for  the  roughness  spectrum. 

In  the  same  figure,  the  received  specular  path  loss,  based 
only  on  the  spherical  spreading,  is  86  dB. 

Figure  2  predicts  a  higher  received  off-specular  level  and  a 
faster  decay  rate  for  a  larger  spectral  power  law.  Over  the 
range  of  arrival  time,  the  selected  wave  numbers  are  in  the  range 
of  about  0.09  to  0.3  (meter  x ) ,  i.e.  magnitudes  of  less  than  one. 
The  selected  spectral  shape  E(k)  =  E  k  m,  indicates  that  the 
case  of  higher  power  law  assumes  a  h?gher  spectral  level  due  to 
the  magnitude  of  k  being  less  than  one.  Therefore,  the  greater 
received  level  for  a  larger  m  value  is  due  to  the  higher  spectral 
level  used.  For  the  same  spectral  level,  E  should  be  adjusted 
for  different  m.  However,  the  higher  decayrate  for  larger  m  is 
physically  caused  by  the  faster  drop  in  the  roughness  spectral 
level  as  the  arrival  time  increases  which  increases  the  magnitude 
of  the  selected  wave  number. 

3 

According  to  measurements  of  Spofford  et  al  ,  the  value  of  m 
is  about  3.  For  this  power  law,  the  off-specular  received  level 
is  found  to  be  independent  of  frequency.  For  m  =  2,  frequency 
dependence  is  depicted  in  Figure  3. 

5 .  Results  for  Sedimentary  Boundary 

Time  spread  catised  by  off-specular  scattering  from  a 
transmitting  rough  boundary  is  investigated  here.  ^For  this  case, 
the  full  scattering  strength  equation  given  by  Kuo-1  is  utilized. 
The  general  equation  is  such  that  an  average  over  two  quadrants 
is  needed.  Otherwise  computations  are  similar  to  the  rigid 
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boundary  case. 

Figure  4  presents  received  levels  from  various  kinds  of 
bottom  boundaries  represented  by  various  combinations  of  bottom 
to  water  density  ratio  /  and  wave  number  ratio  k.  A 
sediment  represented  by  (  f ,  -fe,  )  =  (1.7,  0.8)  was  previously 
used  by  Kuo  because  it  presented  the  best  fit  to  earlier 
experimental  back-scattering  data.  This  sediment  type  together 
with  a  basalt  type  bottom,  gives  received  off-specular  levels 
which  are  very  close  to  those  of  the  rigid  boundary  investigated 
previously.  A  sediment  Soup  type  bottom  (expected  in  the  area 
far  away  from  continental  shelf)  gives  very  low  off-specular 
levels . 

Figure  5  presents  frequency  dependence.  Though  time  decay 
rate  is  higher  in  general,  frequency  dependence  is  similar  to 
that  of  the  rigid  boundary  case. 

Figure  6  depicts  dependence  on  the  grazing  angle.  In  all 
cases,  the  total  transmission  loss  is  geometrically  adjusted  to 
have  the  same  value. 

6 .  Concluding  Remarks 

The  results  presented  here  represent  a  first  attempt  to  model 
the  time  spreading  due  to  scattering  using  perturbation  theory. 
Detailed  comparison  with  experimental  data  is  not  performed. 
However,  private  communication  with  experimentalists  do  indicate 
that  the  predicted  decay  rate  is  representative  for:  some  cases. 
There  are  some  experimental  data  which  show  a  much  slower  decay 
rate.  For  these  cases  further  studies  are  needed. 

The  possible  causes  of  a  much  slower  decay  rate  are  many. 

A  few  possibilities  are  suggested  below. 

1)  Multi-rough  boundaries  associated  with  multi-layers,  can 
cause  multi-off-specular  returns  and  cause  a  slower  decay  rate 
compared  to  a  single  rough  boundary. 

2)  The  combination  of  boundary  and  volumetric  scattering  can 
also  cause  a  slower  decay. 

3)  A  bottom  material  that  supports  shear  can  alter  both  the 
decay  rate  and  the  off-specular  scattering  levels  considerably. 
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